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We determine the continuum limit of the curvature of the pseudocritical line of QCD with n/=2+l 
staggered fermions at nonzero temperature and quark density. We perform Monte Carlo simulations 
at imaginary baryon chemical potentials, adopting the HISQ/tree action discretization, as imple¬ 
mented in the code by the MILC collaboration. Couplings are adjusted so as to move on a line of 
constant physics, as determined in Ref. [I|], with the strange quark mass fixed at its physical 
value and a light-to-strange mass ratio mijma = 1/20. The chemical potential is set at the same 
value for the three quark species, pi = p,s = p. We attempt an extrapolation to the continuum 
using the results on lattices with temporal size up to Lt = 12. Our estimate for the continuum value 
of the curvature k at zero baryon density, k = 0.020(4), is compared with recent lattice results and 
with experimental determinations of the freeze-out curve. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw 


INTRODUCTION 

The phase diagram of strongly interacting matter in 
the temperature (T) - baryon density plane remains a 
challenge for theoretical physics. Although there is little 
doubt that it features a low-temperature hadronic phase, 
with broken chiral symmetry, and a high-temperature de- 
confined phase, with restored chiral symmetry, the ques¬ 
tion about the precise location and the exact nature of 
the transition between these two phases is still open. Yet, 
the answer to this question has many phenomenological 
implications: the region of the phase diagram with high 
T and small baryon density is relevant for the physics of 
the early Universe, whereas the region of low T and high 
baryon density is interesting for the astrophysics of some 
compact objects, but other corners of the phase diagram 
are not less interesting (see Ref. Q for an overview). 

Relativistic heavy-ion collisions provide us with a 
unique opportunity to infer properties of the transition: 
depending on the energy of the ion beams and on the 
mass number of the colliding ions, the fireball generated 
in the collision could fulfill the temperature and baryon 
density conditions under which the deconfined phase ap¬ 
pears as a transient state, before the system freezes 
out into hadrons, which are then detected. Thermal- 
statistical models, assuming approximate chemical equi¬ 
librium at the chemical freeze-out point, are able to de¬ 
scribe the particle yields at a given collision energy in 
terms of two parameters only, the freeze-out tempera¬ 
ture T and the baryon chemical potential pB- The col¬ 
lection of freeze-out parameters extracted from experi¬ 
ments with different collision energy lie on a curve in 


the (T, ^sj-plane, extending up to pB < 800 MeV (see 
Fig. 1 of Ref. [^, or Ref. Q for a recent re-analysis of 
experimental data). 

Quantum ChromoDynamics (QCD) is widely accepted 
as the theory of strong interactions and, as such, must 
encode all the information needed to precisely draw the 
phase diagram in the (T, /rB)-plane. As a matter of fact, 
only some corners of it can be accessed by first-principle 
applications of QCD, in the perturbative or in the non- 
perturbative regime. Here we focus on the lattice ap¬ 
proach of QCD, based on the idea of discretizing the the¬ 
ory on a Euclidean space-time lattice and simulating it 
by Monte Carlo numerical simulations as a statistical sys¬ 
tem, with Boltzmann weight given by exp(—5 'b), where 
Se is the QCD Euclidean action. The region of the phase 
diagram where pbH^T) ^ 1 is within the reach of this 
approach and one can therefore address, at least inside 
this region, the problem of determining the shape taken 
by the QCD pseudocritical line separating the hadronic 
from the deconfined phase. 

There is no a priori argument for the coincidence of 
the QCD pseudocritical line with the chemical freeze- 
out curve: if the deconfined phase is realized in the fire¬ 
ball, in cooling down the system first re-hadronizes, then 
reaches the chemical freeze-out. This implies that the 
freeze-out curve lies below the pseudocritical line in the 
Pb-T plane. It is a common working hypothesis that the 
delay between chemical freeze-out and rehadronization is 
so short that the two curves lie close to each other and 
can therefore be compared. Under the assumptions of 
charge-conjugation invariance at /is =0 and analyticity 
around this point, the QCD pseudocritical line, as well as 
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the freeze-out curve, can be parameterized, at low baryon 
densities, by a lowest-order expansion in the dimension¬ 
less quantity ^J,B as 


Tc{0) 


= 1 — K, 


TilJiB) 


( 1 ) 


where Tc(0) and k are, respectively, the pseudocritical 
temperature and the curvature at vanishing baryon den¬ 
sity. 

Direct Monte Carlo simulations of lattice QCD at 
nonzero baryon density are hindered by the well known 
“sign problem”: Se becomes complex and the Boltzmann 
weight loses its sense. Several ways out of this problem 
have been devised (see Ref. for a review): redesign¬ 
ing the Monte Carlo updating algorithms for a complex 
action reweighting from the ensemble at /tb = 0 Q, 
Taylor expanding the relevant observables around = 0 
and calculating the first coefficients of the series by sim¬ 
ulations at ub = 0 id, using the canonical formu¬ 
lation [H [ij, using the density of states method [l^ 
and simulating the theory at imaginary chemical poten¬ 
tials and performing the analytic continuation to real 
ones [IJSll. 

The numerical evidence gathered so far in QCD with 
Uf = 2-\-l and physical or almost physical quark masses 
points to a scenario with a smooth crossover between 
the hadronic and the deconfined (or chirally symmetric) 
phase at p,B = 0, with a pseudocritical temperature Tc(0) 
of about 155 MeV [E 22|. This crossover behavior should 
persist in some neighborhood of p,B = 0, up to the onset 
of a first-order transition at some value of p,B > 0. 

The state-of-the-art of lattice determinations of the 
curvature k, up to the very recent mpers of Ref. [IS EH, 
is summarized in Fig. 10 of Ref. |19| |: depending on the 
lattice setup and on the observable used to probe the 
transition, the value of k can change even by almost a 
factor of three. The lattice setup dependence stems from 
the kind of adopted discretization, the lattice size, the 
choice of quark masses and chemical potentials, the pro¬ 
cedure to circumvent the sign problem. This dependence 
would totally disappear if, ideally, all groups would use 
the same lattice setup. A contribution to the understand¬ 
ing of the impact of the lattice setjm dependence is pro¬ 
vided in the Appendix B of Ref. [l9|. The dependence 
on the probe observable is, instead, irreducible: since a 
smooth crossover is taking place rather than a true phase 
transition, one cannot define a bona fide order parame¬ 
ter whose behavior would permit to uniquely locate the 
transition point; instead, for any adopted surrogate ob¬ 
servable, a different transition point should be expected, 
at least in principle. 

On the side of the determinations of the freeze-out 
curve, two recent determinations [3, of k, both based 
on the thermal-statistical model, but the latter of them 
including the effect of inelastic collisions after freeze-out. 


give two quite different values of k, each seeming to prefer 
a different subset of lattice results (see Fig. 3 of Ref. [I^ 
for a snapshot of the situation). 

The aim of this work is to contribute to a better under¬ 
standing of the systematics underlying lattice determina¬ 
tions of the curvature k, by corroborating our previous 
determination 0 with an extrapolation to the contin¬ 
uum limit and by comparing it with experimental anal¬ 
yses of the freeze-out curve. 

Our lattice setup is as follows. We simulate the 
HISQ/tree action of the MILC collaboration with 2-|-l 
staggered fermions on lattices with temporal extension 
Lt = 6, 8, 10 and 12 and aspect ratio equal to four. 
We work on the line of constant physics (TCP) as deter¬ 
mined in Ref. [ij, with the strange mass set at the phys¬ 
ical value and the light quark mass fixed at mi = ms/20. 
As discussed in Ref. [H , this amounts to tune the strange 
quark mass until the mass of the fictitious rjss meson 
matches the lowest order perturbation theory estimate 
mrj^g = ^J2m\ — toJ. Consequently within our simula¬ 
tions the pion mass is m-^ ~ 160 MeV. 

We perform simulations at imaginary quark chemical 
potentials, assigning the same value to the three quark 
species, /r; = /Ts = /i, then extrapolate to real chemical 
potentials. Our probe observables are the disconnected 
susceptibility of the light quark chiral condensate and 
its renormalized counterpart. Simulating the theory at 
imaginary chemical potentials poses no restriction on the 
lattice size or in the choice of the couplings. However, the 
periodicity in Im(/r/)/T of the partition function im¬ 
plies that the information gathered outside a narrow in¬ 
terval of imaginary chemical potentials is redundant. For 
the setup with fii = fj,s = n, this interval can be chosen as 
the region 0 < Im(/j,)/T < 7r/3. A safe extrapolation of 
the critical line to real chemical potentials requires that 
it exhibits a smooth dependence on imaginary chemical 
potential over this interval, a condition which must be 
checked to be satisfied a posteriori by our data. 

The preference to the disconnected susceptibility of 
the light quark chiral condensate has multiple motiva¬ 
tions [l|. First of all, for small enough quark masses, its 
contribution to the chiral susceptibility dominates over 
the connected one, which is harder to compute; then, it 
shows a strong sensitivity to the transition; finally, it is 
exempt from additive renormalization, undergoing only 
a multiplicative one. This translate into a very precise 
determination of the critical couplings at imaginary /r, 
which is the main prerequisite of a safe extrapolation to 
the real values of p,. 

There are two main limitations in our setup. The first 
is that we work with a physical strange quark mass, but 
with light quarks a bit heavier than physical ones. Nu¬ 
merical results in nf = 2 indicate a mild dependence of 
the curvature on the quark mass (see the discussion in 
Sec. Ill of Ref. 17|). If the same applies here, as we be¬ 
lieve, our result for k will only slightly underestimate (in 
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absolute value) the true physical curvature. The second 
limitation is that, for the sake of comparison with the 
freeze-out curve, our setup of chemical potentials could 
not be the one which better reproduces the initial condi¬ 
tions of heavy ion collision. In fact, strangeness neutral¬ 
ity would rather impose ^ In general, the setup 
^J-l = ^J-s = fJ'B/S approximates strangeness neutrality at 
low temperatures, while the fig = 0 setup is relevant for 
high enough temperatures. 

It is natural to expect that the effect of taking /ig = /i/ in¬ 
stead of /is = 0 becomes less and less evident when /i //T 
approaches zero, so that the curvature k at zero baryon 
density should not differ too much in the two cases. The 
numerical analysis of Refs. l3“21| has shown that this 
effect is invisible within the accuracy of the lattice setup 
adopted there. 

The paper is organized as follows: in Sec. II we give 
some further details of our numerical simulations. In 
Sec. Ill we show our numerical results for n. Finally, in 
Sec.IV we draw our conclusions. 


SIMULATION DETAILS AND NUMERICAL 
RESULTS 


We perform simulations of lattice QCD with 2-1-1 fla¬ 
vors of rooted staggered quarks at imaginary quark chem¬ 
ical potential. We have made use of the HISQ/tree 
action 23- 2^ a s implemented in the publicly available 
MILC code |27l|. which has been suitably modified by 
us in order to introduce an imaginary quark chemical 
potential /r = /tb/S. That has been done by multiply¬ 
ing all forward and backward temporal links entering the 
discretized Dirac operator by exp(za/j,) and exp(—fa/r), 
respectively: in this way, the fermion determinant is still 
real and positive, so that standard Monte Carlo meth¬ 
ods can be applied. As already remarked above, in the 
present study we have /i = /i; = /ig. This means that 
the Euclidean partition function of the discretized the¬ 
ory reads 


Z = n clet(Djt7,/r])i/\ 

_j „ 


( 2 ) 


q—u,d,s 


where S'gauge is the Symanzik-improved gauge action and 
Dq[U,ij] is the staggered Dirac operator, modified as ex¬ 
plained above for the inclusion of the imaginary quark 
chemical (see Ref. and appendix A of Ref. 261 for 
the precise definition of the gauge action and the covari¬ 
ant derivative for highly improved staggered fermions). 

All simulations make use of the rational hybrid Monte 
Carlo (RHMC) algorithm. The length of each RHMC 
trajectory has been set to I.O in molecular dynamics time 
units. 

We have simulated the theory at finite temperature, 
and for several values of the imaginary quark chemical 


potential, near the transition temperature, adopting lat¬ 
tices of size 16^ X 6, 24^ x 6, 32^ x 8, 40^ x 10 and 48^ x 12. 
We have discarded typically not less than one thousand 
trajectories for each run and have collected from 4k to 
8k trajectories for measurements. 

The pseudocritical point /3c (z^^) has been determined 
as the value for which the renormalized disconnected sus¬ 
ceptibility of the light quark chiral condensate divided by 
exhibits a peak. 

The bare disconnected susceptibility is given by: 

2 

, (3) 


Here Uf = 2 is the number of light flavors and Lg denotes 
the lattice size in the space direction. The renormalized 
chiral susceptibility is defined as: 


Xtren — X^.disc- 


(4) 


The multiplicative renormalization factor can be de¬ 
duced from an analysis of the line of constant physics for 
the light quark masses. More precisely, we have p6| : 


ZM 


miP) 

mi{l3*) ’ 


(5) 


where the renormalization point /3* is chosen such that: 


ri 

a{P*) 


2.37 , 


( 6 ) 


where the function a{l3) is discussed below. In Fig. [T] is 
shown the multiplicative renormalization factor de¬ 
termined in the case when ri is used to set the scale 
(see below). To precisely localize the peak in XJ,ren/T^, a 
Lorentzian fit has been used. For illustrative purposes, in 
Fig. [5] we display our determination of the pseudocritical 
couplings at ^/{ttT) = 0.2/ for all lattices considered in 
this work. The complete collection of results for the dis¬ 
connected susceptibility of the light quark chiral conden¬ 
sate obtained in this work is presented in the Appendix. 

To get the ratios Tc{^)/Tc{Q), we fix the lattice spacing 
through the observables ri and following the discus¬ 
sion in the Appendix B of Ref. (l| . 

For the ri scale the lattice spacing is given in terms of 
the ri parameter as: 


co/(/3) + C2(10//3)/3(/3) 

l + d2(10//3)/2(/3) 


(7) 


with Co = 44.06, C 2 = 272102, d 2 = 4281, ri = 
0.3106(20) fm. 

On the other hand, in the case of the Jk scale we have: 


— 0.05mg 


4 f{l3)+ 4 {10/13) f HP) 
1 + {10 /P)f HP) 


( 8 ) 
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FIG. 1. The multiplicative renormalization factor Zm in the 
case of ri-scale. The renormalization point is /3* = 6.54706. 


with = 7.66, cf = 32911, df = 2388, 0.1738. 

In Eqs. © and ([5]), /(/3) is the two-loop beta function, 

/(/?) = (6o( 10/^))-'^/(2;,^) exp(-/?/(206o)) , (9) 

bo and &i being its universal coefficients. 


Our results are summarized in Table HI For all lattice 
sizes but 24^ x 6 (where we have only one value of /t), 
the behavior of Tc(/r)/Tc(0) can be nicely fitted with a 
linear function in 


Tc{p) 

Tc{0) 


I + Rg 



2 


( 10 ) 


which gives us access to the curvature Rq and, hence, to 
the curvature parameter k = —i?g/(97r^) introduced in 
Eq. (IT|). On the 24^ x 6 lattice the linearity in /j,^ has 
been assumed to hold, in order to extract Rg from the 
only available determination at fi/ijrT) = 0.2i. 

For the sake of the extrapolation to the continuum 
limit, in Fig. [3] we report our determinations of Rq on 
the lattices 24^ x 6, 32^ x 8, 40^ x 10, 48^ x 12, and from 
the two different methods to set the scale, versus l/Tj. 

Within our accuracy, cutoff effects on Rg are negligible, 
so that a constant ht works well over the whole region 
(xl — 0.99), thus including also the smallest 24^ x 6 
lattice. Taking into account the uncertainties due to the 
continuum limit extrapolation. 


K = 0.020(4) . (11) 


Our estimate of the uncertainties for the curvature given 
in Eq. m takes into account both the error in the fit 


FIG. 2. The real part of the renormalized susceptibility of the 
light quark chiral condensate over on the lattices 16® x 6, 
32® X 8, 40® X 10 and 48® x 12 at /r/(7rr) = 0.2i. Full lines 
give the Lorentzian fits near the peaks. The temperature has 
been determined from the ri scale. 


minimization and the choice of the minimization func¬ 
tion. We stress, however, that if we exclude from the fit 
the value on the lattice with the smallest Lj, i.e. the 
rightmost points in Fig. [3l the extrapolation to the con¬ 
tinuum becomes largely undetermined. Indeed, with the 
values of Rg obtained in the present work (see Table HIl) . 
the fit with a constant is rather stable even if Lj = 6 is 
excluded, but the fit with a linear function in 1 /L^ in the 
latter case gives a much smaller value of the curvature k, 
though with a large uncertainty (see Table IIII|) . 


CONCLUSIONS AND DISCUSSION 

We have studied QCD with Uf = 2 + 1 flavors dis¬ 
cretized in the HISQ/tree rooted staggered fermion for¬ 
mulation and in the presence of an imaginary baryon 
chemical potential, with a physical strange quark mass 
and a light-to-strange mass ratio mi/rris = 1/20, and 

fl = 

We have estimated, by the method of analytic con¬ 
tinuation, the continuum limit of the curvature of the 
pseudocritical line in the temperature - baryon chemical 
potential, defined in Eq. O- The observable adopted 
to identify, for each fixed p,, the crossover temperature 
has been the disconnected part of the renormalized sus¬ 
ceptibility of the light quark chiral condensate, in units 
of the squared temperature. This observable is conve¬ 
nient for many reasons: it dominates, for small enough 
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TABLE I. Summary of the values of the ratio Tc{fJ-)/Tc{0) for 
the imaginary quark chemical potentials fj, considered in this 
work. The data for /r = 0 on the 24® x 6, 32® x 68 and 48® x 
12 lattices have been estimated from the disconnected chiral 
susceptibilities reported respectively on Tables X, XI and XII 
of Ref. P]. The datum for fj, — 0 on the 40® x 10 lattice has 
been estimated from the disconnected chiral susceptibilities 
reported on Table XI of Ref. [^. The values of /Tc{0) 
evaluated fixing the lattice scale by ri and /k are reported, 
respectively, in the third and in the fourth column of the table. 


lattice 

^J./{nT) 

T 0 M/T 4 O) 

(n scale) 

T4b)/To(0) 
(/k scale) 

16® X 6 

O.lSi 

1.038(13) 

1.043(14) 


0 . 2 i 

1.063(15) 

1.070(15) 


0.25i 

1.085(16) 

1.095(18) 

24® X 6 

0 . 2 i 

1.061(9) 

1.067(10) 

32® X 8 

0.15i 

1.054(7) 

1.059(8) 


0 . 2 i 

1.066(10) 

1.071(11) 


0.25i 

1.117(10) 

1.126(10) 

40® X 10 

0.15i 

1.023(23) 

1.024(24) 


0 . 2 i 

1.075(14) 

1.079(15) 


0.25i 

1.102(15) 

1.107(15) 

48® X 12 

0.15i 

1.013(31) 

1.013(33) 


0 . 20 i 

1.051(14) 

1.052(15) 


0.25i 

1.094(26) 

1.097(25) 


quark masses, the whole light chiral susceptibility, which 
would be much harder to implement; it undergoes only 
a multiplicative renormalization; it is strongly sensitive 
to the transition, thus allowing precise determinations of 
the pseudocritical temperatures. 

We have found that, within the accuracy of our deter¬ 
minations, cutoff effects on the curvature are negligible 
already on the lattice with temporal size Lt = 6. Our 
determination of the curvature parameter, k= 0.020(4), 
is indeed compatible with the value quoted in our previ¬ 
ous paper [l^, k= 0.018(4), without the extrapolation to 
the continuum. 

It is interesting to extrapolate the critical line as de¬ 
termined in this work to the region of real baryon density 
and compare it with the freeze-out curves resulting from 
a few phenomenological analyses of relativistic heavy-ion 
collisions. This is done in Fig.|4j where we report two dif¬ 
ferent estimates. The first is from the analysis of Ref. Q, 
based on the standard statistical hadronization model, 
where the freeze-out curve is parametrized as 

Tcifis) =a-bn\- cfi% , ( 12 ) 

with a = 0.166(2) GeV, b = 0.139(16) GeV\ and 
c = 0.053(21) GeV~^. The second estimate is from 
Ref. [ 2 ^ and is based on the analysis of susceptibilities of 
the (conserved) baryon and electric charges. In fact, our 
critical line is in nice agreement with all the freeze-out 
points of Refs. im. In particular, using our estimate 
of the curvature, Eq. (|TT1) . we get b = 0.128(25) GeV“^, 
in very good agreement with the quoted phenomenolog¬ 
ical value. The significance of the comparison presented 



FIG. 3. Determinations of the curvature Rq on the lattices 
24® X 6, 32® X 8, 40® x 10, 48® x 12, and from the two different 
methods to set the scale, versus l/L?. Data points related 
with the /k scale setting have been slightly shifted along the 
horizontal axis for better readability. The dashed horizontal 
line gives the result of the fit to all data with a constant; the 
solid horizontal lines indicate the uncertainty on this constant 
(95% confidence level). 

in Fig. |4l with special reference to the question whether 
the pseudocritical line lies indeed above the freeze-out 
curve, can be increased at the (nonnegligible) price of 
reducing the uncertainties on T(,(0) and on k. 

Some caveats are in order here. We do not expect our 
critical line to be reliable too far from fi = 0: as a rule 
of thumb, we can trust it up to real quark chemical po¬ 
tentials of the same order of the modulus of the largest 
imaginary chemical potential included in the fit ^■e. 
|^|/(7rT) = 0.25. This translates to real baryon chem¬ 
ical potentials in the region ^ 0.4 GeV. Moreover, 
the effect of taking /is = fii instead of /is < /i; should 
become visible on the shape of the critical line as we 
move away from /i = 0 in the region of real baryon den¬ 
sities, thus reducing further the region of reliability of 
our critical line. So, from a prudential point of view, the 
agreement shown in Fig. 3] could be considered the fortu¬ 
nate combination of different kinds of systematic effects. 
We cannot however exclude the possibility that the mes¬ 
sage from Fig.|4]is to be interpreted in positive sense, i.e. 
the setup we adopted and the observable we considered 
may catch better some features of the crossover transi¬ 
tion, thus explaining the nice comparison with freeze-out 
data. Indeed, our result for the continuum extrapolation 
of the curvature k is in fair agreements with the recent es- 
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TABLE II. Summary of determinations of the curvature Rq for all values of Lt considered in this work, and from the two 
different methods to set the scale. 


Lt 

6 

8 

10 

12 

Rq (ri scale) 

K (ri scale) 

-1.466(306) 

0.017(3) 

-1.902(192) 

0.021(2) 

-1.685(294) 

0.019(3) 

-1.337(410) 

0.015(5) 

Rq {/k scale) 

K {rK scale) 

-1.646(336) 

0.019(4) 

-2.041(206) 

0.023(2) 

-1.769(309) 

0.020(3) 

-1.394(415) 

0.016(5) 


TABLE III. Summary of the fit of the curvature k with the function K{Lt) = k + Aj, with A taken equal to zero or left 
free. The first column specifies the values of Lt included in the fit, the last column the reduced ■ The uncertainties on the 
fit parameters are obtained with 70% confidence level. 


Lt included 

K 

A 

xi 

6, 8, 10, 12 

0.01991(114) 

0. 

0.76 

6, 8, 10, 12 

0.02014(449) 

-0.015(259) 

0.89 

8, 10, 12 

0.02048(127) 

0. 

0.80 

8, 10, 12 

0.01182(748) 

0.669(559) 

0.14 



FIG. 4. Tc{^j.b) versus (units in GeV). Experimental 
values of TcIhb) are taken from Fig. 1 of Ref. (black 
circles) and from Fig. 3 of Ref. (green triangles), for 
the standard hadronization model and for the susceptibili¬ 
ties of conserved charges respectively. The dashed line is a 
parametrization corresponding to Tc{^jlb) = T'c(O) — bfi% with 
rc(0) = 0.154(9) GeV and b = 0.128(25) GeV^ . The solid 
lines represent the corresponding error band. 


timates in Ref. [20j , where both setup fig = Mi S'lid Ms = 0 
were adopted, and Ref. [21|, where the strangeness neu¬ 
tral trajectories were determined from lattice simulations 
by imposing (ns) = 0. 
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Appendix: Summary of data for the disconnected 
susceptibility of the light quark chiral condensate 

In this Appendix (Tables IIVI through IXVI|) , we sum¬ 
marize all the results for the disconnected chiral suscep¬ 
tibility obtained in our simulations for each considered 
value of the coupling /3 and for the corresponding phys¬ 
ical temperature, as determined from the two different 
procedures to set the scale. 
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TABLE IV. Data for the disconnected chiral susceptibility on the 16® x 6 lattice at ^/{'kT) = O.lSOi: the second column gives 
Xi.disc defined in Eq. ([3|, the fourth and sixth columns give Xi.ren/T^, with Xi,ren defined in Eq. (|4]) and the temperature T 
(columns three and five, respectively) determined by the two different methods to set the scale. 


p 


T(MeV) 

(ri scale) 

Xi.ren/T^ 

r(MeV) 
(/k scale) 

Xi.ren/T^ 
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0.939 (0.114) 

147.3 

25.5 (3.1) 
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27.7 (3.4) 

6.050 
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166.7 
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6.195 

1.093 (0.107) 

178.1 

34.1 (3.3) 

170.2 

35.6 (3.5) 

6.200 

0.944 (0.100) 

178.9 

29.5 (3.1) 

171.1 

30.8 (3.3) 


TABLE V. Data for the disconnected chiral susceptibility on the 16® x 6 lattice at /i/(7rT) = 0.200i (legend as in Table HVTl . 
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TABLE VI. Data for the disconnected chiral susceptibility on the 16® x 6 lattice at ^/{tvT) = 0.250i (legend as in Table HVl) . 
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TABLE VIII. Data for the disconnected chiral susceptibility on the 32® x 8 lattice at ^/(tvT) = O.lSOi (legend as in Table ITVT) . 


p 

X^idisc 

T(MeV) 

(ri scale) 

Xi,ren/T^ 

r(MeV) 
(/k scale) 

Xi.ren/T^ 

6.415 

0.646 (0.073) 

165.5 

39.7 (4.5) 

160.9 

40.1 (4.6) 

6.424 

0.730 (0.047) 

167.0 
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55.4 (4.4) 
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0.776 (0.060) 

174.6 

48.6 (3.8) 
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48.8 (3.8) 


TABLE IX. Data for the disconnected chiral susceptibility on the 32® x 8 lattice at /r/(7rT) = 0.200i (legend as in Table HVl) . 
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TABLE X. Data for the disconnected chiral susceptibility on the 32® x 8 lattice at ^/(ttT) = 0.250* (legend as in Table HVll . 


p 

X^jdisc 

T(MeV) 

(ri scale) 

Xi,ren/T^ 

r(MeV) 
(/k scale) 

Xi.ren/T^ 

6.420 

0.613 (0.055) 

166.4 

37.7 (3.4) 

161.7 

38.2 (3.4) 

6.455 

0.555 (0.044) 

172.1 

34.6 (2.7) 

167.7 

34.8 (2.7) 

6.490 

0.792 (0.053) 

178.1 

49.9 (3.4) 

173.8 

50.0 (3.4) 

6.525 

0.753 (0.051) 

184.2 

47.8 (3.3) 

180.2 

47.8 (3.3) 

6.560 

0.568 (0.030) 

190.6 

36.5 (2.0) 

186.7 

36.3 (1.9) 


TABLE XL Data for the disconnected chiral susceptibility on the 40® x 10 lattice at ^/(tvT) = 0.150* (legend as in Table ITVT) . 


P 


T(MeV) 

(ri scale) 

Xi.ren/T^ 

r(MeV) 
(/k scale) 

X*.ren/r^ 

6.550 

0.444 (0.043) 

151.0 

44.4 (4.3) 

147.9 

44.3 (4.3) 

6.606 

0.570 (0.070) 

159.4 

57.9 (7.1) 

156.5 

57.4 (7.0) 

6.648 

0.524 (0.037) 

165.9 

53.6 (3.8) 

163.2 

53.0 (3.8) 

6.690 

0.525 (0.048) 

172.7 

54.1 (4.9) 

170.2 

53.3 (4.8) 

6.732 

0.402 (0.035) 

179.8 

41.7 (3.6) 

177.4 

41.0 (3.5) 


TABLE XII. Data for the disconnected chiral susceptibility on the 40® x 10 lattice at ^,/{tvT) — 0.200* (legend as in Table HVll . 


P 

X^jdisc 

T(MeV) 

(ri scale) 

Xi,ren/T^ 

r(MeV) 
(/if scale) 

Xi.ren/T^ 

6.575 

0.420 (0.042) 

154.7 

42.3 (4.3) 

151.7 

42.0 (4.2) 

6.600 

0.425 (0.047) 

158.4 

43.1 (4.8) 

155.5 

42.7 (4.7) 

6.630 

0.407 (0.043) 

163.1 

41.6 (4.4) 

160.3 

41.1 (4.3) 

6.655 

0.525 (0.057) 

167.0 

53.8 (5.8) 

164.4 

53.1 (5.8) 

6.680 

0.550 (0.067) 

171.1 

56.7 (6.9) 

168.5 

55.8 (6.8) 

6.730 

0.481 (0.051) 

179.4 

49.9 (5.3) 

177.1 

49.0 (5.2) 

6.775 

0.297 (0.037) 

187.3 

31.0 (3.9) 

185.1 

30.4 (3.8) 


TABLE XIII. Data for the disconnected chiral susceptibility on the 40® x 10 lattice at fx/{'KT) = 0.250* (legend as in Table ITVTl . 


P 

X^^disc 

T(MeV) 

(ri scale) 

Xi.ren/T^ 

r(MeV) 
(/if scale) 

Xi.ren/T^ 

6.645 

0.386 (0.060) 

165.4 

39.5 (6.1) 

162.7 

39.0 (6.1) 

6.664 

0.368 (0.041) 

168.5 

37.8 (4.2) 

165.8 

37.3 (4.1) 

6.714 

0.552 (0.064) 

176.7 

57.2 (6.6) 

174.3 

56.3 (6.5) 

6.764 

0.421 (0.059) 

185.3 

43.9 (6.1) 

183.1 

43.0 (6.0) 

6.810 

0.193 (0.051) 

193.6 

20.2 (5.3) 

191.5 

19.8 (5.2) 


TABLE XIV. Data for the disconnected chiral susceptibility on the 48® x 12 lattice at ^/{'kT) — 0.150* (legend as in Table ITVll . 


P 

X^ jdisc 

r(MeV) 

(ri scale) 

Xi,ren/T^ 

r(MeV) 

(/if scale) 

Xi,ren/T^ 

6.700 

0.242 (0.052) 

145.3 

36.1 (7.7) 

145.3 

36.1 (7.7) 

6.730 

0.374 (0.047) 

149.5 

55.9 (7.1) 

149.5 

55.9 (7.1) 

6.780 

0.433 (0.063) 

156.8 

65.3 (9.5) 

156.8 

65.3 (9.5) 

6.815 

0.475 (0.111) 

162.1 

71.9 (16.8) 

162.1 

71.9 (16.8) 

6.830 

0.382 (0.036) 

164.4 

57.9 (5.5) 

164.4 

57.9 (5.5) 

6.880 

0.466 (0.068) 

172.3 

71.1 (10.4) 

172.3 

71.1 (10.4) 

6.930 

0.247 (0.027) 

180.5 

37.8 (4.1) 

180.5 

37.8 (4.1) 

6.980 

0.191 (0.032) 

189.1 

29.3 (4.9) 

189.1 

29.3 (4.9) 
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TABLE XV. Data for the disconnected chiral susceptibility on the 48^ x 12 lattice at /r/(7rT) = 0.200i (legend as in Table ITVTl . 


p 

X^>disc 

T(MeV) 

(ri scale) 

Xi,ren/T^ 

r(MeV) 
(/if scale) 

Xi.ren/T^ 

6.770 

0.376 (0.049) 

155.3 

56.5 (7.4) 

153.5 

55.3 (7.3) 

6.804 

0.357 (0.033) 

160.4 

54.0 (5.0) 

158.6 

52.8 (4.8) 

6.864 

0.487 (0.054) 

169.7 

74.2 (8.2) 

168.1 

72.2 (8.0) 

6.924 

0.249 (0.028) 

179.5 

38.1 (4.3) 

178.1 

37.0 (4.2) 

6.975 

0.129 (0.034) 

188.3 

19.8 (5.2) 

186.9 

19.2 (5.1) 


TABLE XVI. Data for the disconnected chiral susceptibility on the 48® x 12 lattice at fi/ijiT) = 0.250i (legend as in Table HVTl . 


P 


T(MeV) 

(ri scale) 

Xi.ren/T^ 

r(MeV) 
(/if scale) 

X*.ren/T^ 

6.650 

0.274 (0.044) 

138.5 

40.4 (6.5) 

136.3 

39.9 (6.4) 

6.750 

0.175 (0.054) 

152.4 

26.3 (8.1) 

150.5 

25.8 (8.0) 

6.864 

0.300 (0.042) 

169.7 

45.7 (6.3) 

168.1 

44.5 (6.2) 

6.900 

0.418 (0.043) 

175.6 

64.0 (6.6) 

174.1 

62.2 (6.4) 

6.950 

0.278 (0.037) 

183.9 

42.6 (5.6) 

182.6 

41.4 (5.4) 

6.975 

0.224 (0.025) 

188.3 

34.4 (3.8) 

186.9 

33.4 (3.7) 

7.025 

0.153 (0.018) 

197.2 

23.6 (2.7) 

196.0 

22.8 (2.6) 












